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We report a two-parametric irreducible infinitely dimensional representation of the Lax integrability condi¬ 
tion for the fermi Hubbard chain. Besides being of fundamental interest, hinting on possible novel quantum 
symmetry of the model, our construction allows for an explicit representation of an exact steady state many- 
body density operator for non-equilibrium boundary-driven Hubbard chain with arbitrary (asymmetric) particle 
source/sink rates at the letf/right end of the chain and with arbitrary boundary values of chemical potentials. 
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Introduction.- The one-dimensional (ID) fermi Hubbard 
model is one of the key_paradigms of exactly solvable models 
in theoretical physics [l|]. Exact integrability is a remarkable 
exceptional property of a set of classical and quantum mod¬ 
els, playing a prominent role in theoretical physics. Integra¬ 
bility is essentially the only assumption-free handle on non- 
perturbative physics of strong interactions. Exact solution of 
2D Ising model formed a basis for the theory of phase transi¬ 
tions, solution of the Korteweg de Vries equation gave rise to 
soliton theory, solution of the XYZ spin chain model showed 
a breakup of the scaling relations, and solution of ID Hubbard 
model laid a milestone in understanding the problem of itin¬ 
erant electrons in strongly correlated systems. The integrabil¬ 
ity comprises an existence of a Yang Baxter structure, or Lax 
structure, consisting of R-matrix and Lax operators, depend¬ 
ing on a free spectral parameter. The Lax operator generally 
provides a crucial ingredient of integrability, an infinite set of 
mutually commuting local integrals of motion, through e.g. a 
construction of the commuting parameter-dependent transfer- 
matrices M- The one-dimensional Hubbard model stays 
aside in the list of integrable systems, as, since its discovery, 
it continues to raise many deep fundamental questions. Two 
decades after the discovery of the nested Bethe Ansatz for 
Hubbard model SS a Yang-Baxter formulation of the prob¬ 
lem has been found by Shastry |0, H]. Remarkably, Shastry’s 
R-matrix has several unique properties EJ not shared among 
most of other known integrable systems. 

A recent progress has allowed to push the concepts of in¬ 
tegrability of quantum systems beyond the thermal equilib¬ 
rium: exact nonequilibrium steady states (NESS) for a set 
of quantum chains driven far from equilibrium by local in¬ 
coherent markovian noise, were calculated US The non¬ 
equilibrium integrability, it turned out, results from an ad¬ 
ditional degree of freedom, hidden in the standard Yang- 
Baxter structure for its ‘equilibrium’ counterpart. This addi¬ 
tional degree of freedom is a representation parameter (be¬ 
ing a complex number in non-equilibrium setting and con¬ 
nected with the dissipative coupling at the model’s bound¬ 
aries), gives rise to non-unitary representations of quantum 
symmetries of these models and allows to formulate a Ma¬ 
trix Product Ansatz for a NESS. Moreover, these representa¬ 


tions are infinite-dimensional, unlike their ‘equilibrium’ coun¬ 
terparts. 


Existence of an extra representation parameter in the Yang- 
Baxter construction (alongside with the inherited spectral pa¬ 
rameter) not only allowed to find exact non-equilibrium gen¬ 
eralizations of well-known integrable models, but gave rise to 
an intrinsically new infinite set of conservation laws for their 
equilibrium counterparts, through a two-parameter commut¬ 
ing family of operators, generated by monodromy matrix el¬ 
ements. Application of these new conserved quantities has 
allowed to solve a long-standing problem of high-temperature 
ballistic transport in the gapless anisotropic Heisenberg spin- 
1/2 ( XXZ ) chain by establishing a rigorous lower bound on 
the corresponding spin Drude weight SEMI- However, 
for the Hubbard model such a Lax representation was lack¬ 
ing. 


Here we report a Lax structure for the Hubbard model, con¬ 
taining both, a free spectral parameter and an extra (com¬ 
plex) representation parameter. As in previously solved ex¬ 
amples {XXZ model), it is infinite-dimensional and gives rise 
to a two-parameter commuting family of operators generated 
by a monodromy matrix, yielding a new set of conserva¬ 
tion laws. For a fixed value of the spectral parameter, the 
new Lax structure yields an appealing factorization of our 
previous result Q. Remarkably, unlike in the previously 
solved XXZ model, the presented non-equilibrium Lax op¬ 
erator is different from the known one Q] and cannot be re¬ 
duced to the latter by means of auxiliary space truncation and 
gauge transformations. Thus it is not a generalization (like 
in other integrable example El where the non-equilibrium 
Lax operator can be encoded within the so-called universal R- 
matrix of the quantum group U q { sl(2)) symmetry of the XXZ 
model), and it reveals a fundamentally different symmetry of 
the ID Hubbard model. The new Lax solution yields a spec¬ 
trum of potential applications, from quantum transport in con¬ 
densed matter to AdS/CFT duality and integrable N - 4 su¬ 
per Yang-Mills theory lUt l. In the present communication we 
show how our generalised Lax operator can be employed to 
solve asymmetrically boundary driven Hubbard chains with 
arbitrary source/sink rates and boundary chemical potentials, 
where the non-vanishing value of the spectral parameter cor- 
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responds to the asymmetry of the driving. 

We consider a fermi Hubbard chain on n sites formulated 
in terms of a spin 1 /2 ladder. Let TL = C 2 ® C 2 be a lo¬ 
cal physical Hilbert space supporting two independent sets of 
Pauli matrices ct s ,t', s,t e ST '■= {+, —,0, z}, cr° — r° = 1. 
Embedding 7/ p into the Hilbert space of an n-site ladder 7/®, 
one defines local spin operators cr s -, r'-, for j e {1,..., n }. The 
Hubbard hamiltonian with arbitrary boundary chemical poten¬ 
tials yUL/R 112 111 then reads 


n-l 

H = ^ hjj +1 + /;l + Iir, 
j= i 

: = /! U + h h + (°1 T 1 + °1 t 2) - 

^L/R := \o\i&in + ^ (°1/« + ^/a) • 


( 1 ) 

( 2 ) 

(3) 


where h^ 2cr|<x 2 + 2 h\ 2 := 2r|r, + 2 t~t;( are the 
free hoping operators of the respective particle species. Note 
that (1,2) designates a generic pair of neighboring sites. Di¬ 
mensionless interaction parameter u = U /(2f h ) contains stan¬ 
dard Hubbard interaction U and hopping amplitude f h . The 
standard fermionic Hamiltonian H — - .c s j+l +H.c.) + 

2mZ ; («t, 7 - 5)(«u “ 2 ) + + »i,t - 1) + Pr(mt,« + 

W X,« _ 1) is reconstructed via Jordan-Wigner transformation 

C TJ = ^ 7-107 and n,j := cl jC sJ , where 

P { p := (Tj • • • cj, Pp := Tj • • • Ty. Let us define the spin-flip 
(or particle-hole) operator G, i.e. permutation between cr and 
r spins (or fermion species), as Gcr s G = r v , G 2 = 1. Clearly, 
Gh‘{ 2 G = h] v and G/i u G = /q, 2 , G//G = H. 

The Lax operator- We introduce an infinitely dimensional 
auxiliary Hilbert space, 74 = LinearSpan{|/?), p e W), whose 
orthonormal basis is conveniently labelled by vertices of a 
graph (see Fig. Q). = {0 + , ± + , ±“, l - , 1 + , | + , §“, 2", 2 + ...). 

We extend the definition of the spin-flip G over 7/ a as a diag¬ 
onal reflection of the graph, Glfc*) = Ik*), G|k+~ ) = |k+2 ), 
k e Z + . Here and below we shall use bold-roman letters to 
designate operators which are non-scalar over 7/ a . We begin 
our analysis with a simple observation: 

Lemma: Assume there exist operators S, S,S, T, T, T e 
End(7/ a ® 7/p), and X, Y 6 End(7/ a ) (acting as scalars over 
7/p), satisfying 


[/^ 2 ,SiXS 2 ] = SiXS 2 - S,XS 2 , (4) 

[^ 2 ,T,XT 2 ] =T 1 XT 2 -TiXt 2 , (5) 

ST + TS - ST - TS = [Y - wcrV, ST], (6) 

[S, T] = 0, (7) 

[X, Y] = 0. (8) 


Subscripts, like in S ; , indicate independent local physical 
spaces pertaining to sites j. Then, one can define a Lax oper¬ 
ator and its ‘derivative’ L. L e End(7/ a ® 7/ p ) as 


L = STX, (9) 

L = i(ST + TS + ST + TS-{Y,ST})X, (10) 


such that a so-called Sutherland-Shastry (or generalized [7j, 
22;] local operator divergence [230 gLOD) condition holds 

[A u ,L 1 L 2 ] = (L 1 +YL 1 )L 2 -L 1 (L 2 +L 2 Y). (11) 


The proof is a straightforward insertion of (19110b into Eq. (fTTb 
followed by subsequent application of identities ®® observ¬ 
ing the definition ®. 

We continue by deriving an explicit closed form representa¬ 
tion of algebraic identities ®®. Assuming the spin-flip sym¬ 
metry 

GSG = T, GSG = T, GSG = T, [G,X] = [G, Y] = 0, 

( 12 ) 

and writing out the components S = Yjs^j S'cr 1 , T = 
TV, and similarly for S,S, T, T, we find that Eqs. ® 
and ® are equivalent, Eq. ([6]) is invariant under G, while 
Eq.® implies [S',T'] = 0. 

Eqs. m are in fact just a particularly ‘decorated’ 6-vertex 
Yang-Baxter equations for free fermion (or XX) chains. We 
shall thus make an ansatz for S', T' in which each square pla- 
quette [k + ,k + j ,k + \ ,k+\~] of the graph spans a pair of 
representations of a free fermion algebra (see Fig®, namely 
requiring that {S + , S } (and similarly for {T + ,T - } via (fl2b ) is 
in the center of the algebra generated by S s , T' 

[{S + ,S },S 9 ] = [{S + ,S~),T'] =0, s,tej. (13) 


One finds that these conditions are fulfilled by an ansatz 

OO 

s + = V 2 ^(r><k+i + i + ik+i + ><k+ri), G4) 

k= 0 

OO 

s- = V2^(-i/(|k+r><k + | + |k+r><k+i + |), 

k= 0 

OO 

S° = YjP k+ )< 2k+ \ + |2k+i + )<2k+I + | 

k=0 

+\2k+r)(2k+r\ + \2k+p)(2k+p\) 

OO 

+AjP{\2k-p)(2k- k+ \ + \2k-)(2k -\), 

k= 1 

OO 

s z = ^(|2k-l + )<2k-l + | + \2k-p)(2k-p\ 

k=l 

+\2k~)(2k~\ + \2k+p)(2k + p\) 

OO 

+4^] (\2k+p)(2k+p\ + \2k+\-)(2k+\-\), 

k =0 

where A e C is a free parameter. Eqs. (Il2l) imply definition of 
another set of auxilliary fermi operators T' = GS'G, such that 
Eq.® is satisfied. 

Further, one finds a remarkably consistent ansatz for the 
‘interaction’ operator X coupling the neighboring plaquettes: 

OO 

x = |0 + ><0 + | + £(-l)* 2 \k v )xf{k v '\ 

k—l V. !/£(-,+) 
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+<o2(-i) i (i*+r><*+ri+i*+n<*+ri). w 

k =0 

where X* - {Xj' v } VjV 'g{—,+) are still unknown 2x2 matrices 
and weCis another free parameter. Namely, Eq. (0J) yields 
a system of linear equations for auxiliary operators S S X, XS', 
with a unique solution parametrised by X( : , to, A: 

oo 

S + X = -2 V 2 )<k+ i + |, (16) 

k= 1 

oo 

S-X = -2V2j]x,- + |k + )<k-i[, 

Ai=l 

OO 

XS + = 2 V2^](-l)% + -|k-^><k + | 

fc=l 

OO 

XST = -2 V2 ^ Xj~ + \k+j + )(k~\, 

k= 1 

00 

S°X = XS° = 2 ^](w|2k- l + )(2k-1 + | - m|2r)<2r| 

£=1 

-Xi + _ I |2k-i + )<2k-i + |-X 2 7|2k-i-)<2k-i“|) 

OO 

+2d^(-w|2k + )<2k + |+X 2 7 +1 |2k+I')<2k+ri), 

k=0 

OO 

S Z X = XS Z = 2^(w|2k+r><2k+l“| - w|2k + )(2k + | 

k= 0 

+X 2 7|2k+i + )<2k+i + |+X 2 7 +1 |2k+i-)<2k+i _ |) 

OO 

+24 ^](«|2k- l + )(2k-1 + | - X 2 7|2k-0<2k-1"|). 

*=1 

Assuming X to be invertible (i.e., to 4 0, det X\, 4 0) and 
plugging expressions (fl6l) to the remaining identity <(6ji result 
in (i) a unique consistent expression for the ‘spectral’ operator 

Y 


Y = -24n^](|k + ><k + | + |fc+r><k+r|), (17) 

k=0 

which clearly commutes with X, as required by ©, and (ii) 
recurrence relations for the matrix elements of X*: 77 = 
X~~ - uto, = X* + - uto( 1 - A 2 ), and detX* = -to 2 , while 
also fixing the initial condition Xq + = 1, X~" - -to 2 , yielding 

vti \ (~(u + ku)to 1 - (to + ku)co(l - d 2 )\ 

W0J) = 1 -kuto 1 - kuu>{\ - A 2 ) I' (18) 

Note that X7 + /Xt~ can be chosen freely exploring a gauge 
freedom Ik*) —> 7 1 l^ ± )> k = 1,2... We have thus con¬ 
structed two-parameter representation of the Lax operator 
L(A, to) = S(4)T(4)X(4, to) satisfying gLOD We pro¬ 
pose to call A a spectral parameter and to a representation 
parameter. Remarkably, our representation is genetically of 
infinite-dimension, for any nonzero u, and can only be trun¬ 
cated to 4k-dim. span of first k plaquette states at special 



FIG. 1. (Color online) Diagrammatic representation of factors of the 
Lax operator where auxiliary states are labelled by vertices < V. Di¬ 
agrams for T 1 are obtained by reflection of those of S‘ across the 
diagonal. Red/blue arrows indicate offdiagonal transitions with am¬ 
plitude ± V2. Red, blue, green, black, open points represent diagonal 
multiplications by to,-to, oc A, 1,0, respectively, and brown circles 
represent multiplications by 2 x 2 matrices Xj. 


points along algebraic curves Xf + (A, to) = 0 in the A - to 
plane. On the other hand, we have checked that the Shastry’s 
Lax operator |0], together with appropriate auxiliary opera¬ 
tors, form a 4-dim. representation of the algebra (|4][8]> as well. 
So one can be tempted to think that our Lax structure, for 
X^ + (A, to(A)) - 0, and Shastry’s Lax structure can be equiva¬ 
lent; however, so far, we could not find any correspondence. 

Lax form of NESS for asymmetric boundary driving.- As 
an application of the novel Lax operator we consider a marko- 
vian master equation dp,/dr = lip, for an open Hubbard chain 
with Hamiltonian <[T]) and driven by pure source/sink at the 
left/right ends with the non-negative rates I L / R : 


Ip = -i [H,p] + (r L (iv r + i) T ;) + r R (i7- + D T -j) P , ( 19 ) 


where Dip = 2 Lplf - {L 1 L,p) is a Lindblad dissipator 124 , 
2511 . a linear map over Endl'K*"). We shall now prove the 
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following explicit result on NESS poo'- 

Theorem: The unique 0 fixed point £poo — 0 of boundary 
driven Hubbard chain reads 

Poo = (tiR) l R, = (20) 

where O = Q(A, to) is a highest-weiglit transfer matrix 

Q = <0 + |Li(d, lo)L 2 (A, to)--- LJA, w)|0 + ) (21) 


and M is a diagonal operator 

M = MiM 2 ---M„, Mj = exp (rjicr) + rp) (22) 

with r) — j logE L /r R and parameters A,to are related to co¬ 
herent and incoherent biases 


Tl ~ Tr - i(pL + Pr) 
Tl + Tr - i(pL - Pr) 


w - - (pl - Pr + i (Tl + Tr)) . 


(23) 

Proof. Let us now invoke two copies of the auxiliary space 
and define operators S, T, S', T' e End(77 a ® 7Y a ® 77 p ) as 


S = 2s i ®l a ®<r , > 

S 

S' = 2 la 9 S ' ® ( tr s ) T , 


T = ^ T' ® l a ® r', and 

t 

T' = Yj la®T' ®(r') 7 '. 


() T denotes the matrix transposition and S the complex con¬ 
jugation, i.e. replacement A, to, —> A, to, and similarly for 
S, S, S', S', T, T, T', T', and X,X',Y,Y' e End(77 a ® 77 a ). In 
fact, the primed operators S',S',S',T',T',T',X', Y' gener¬ 
ate a conjugate representation of the algebra (|4][8j. Noting 
[h\ 2 ,M\M 2 ] = 0 and the Jacobi identity one finds that the 
following double auxiliary operators 


L j = LjL'jMj, Lj = (LjVj - L jL'j)Mj, Y = Y - Y', (24) 
also respect gLOD ( ITTI) . resulting in the telescoping series 

n— 1 

^ L 1 L 2 • • ■ L„j = (ILi + {Y, Li})L2 • • ■ L„ 

j= 1 

- L] • • -l^ n -\0L n + {Y,lL n }). (25) 


Double Lax operator expresses NESS in a compact form 


R = <0 + ,0 + |L 1 L 2 ---L„|0 + ,0 + ), (26) 

hence the fixed point condition JjR = 0 becomes, after apply¬ 
ing (l25l) to [//. R\, equivalent to a pair of equations for ultralo¬ 
cal operators at the boundary physical sites 

<0 + , 0 + | (iELliV + i) T+ )L + L + LY + [h L , L]) = 0, 

(iE R (iV + B t -)L - L - YL + [/ iR , L]) |0 + , 0 + ) = 0, (27) 

where boundary interactions with fields, Pl/r, are defined in 
©. Using explicit forms (1141118b and in particular X|0 + , 0 + ) = 
X'|0 + ,0 + ) = |0 + ,0 + ), each of Eqs. (|27T> results in dimTfp X 
dim 7Y p = 16 equations for (bra/ket) vectors from 77 ;l ® 7f a , 


most of them trivially satisfied, whereas the non-trivial ones 
being equivalent to conditions (l23l . 

Discussion.- We have derived an infinitely dimensional 
irreducible representation of Lax operator and Sutherland- 
Shastry compatibility condition and shown how it can be 
employed to yield exact NESS of asymmetrically boundary 
driven Hubbard chain with arbitrary boundary chemical po¬ 
tentials. One is thus able to explicitly compute physical ob¬ 
servables in NESS. For example, linear dependence of the 
amplitudes (fTSt on auxiliary state k immediately yields, sim¬ 
ilarly as in Heisenberg chain 10, a universal scaling of the 
spin/charge currents J ~ n 2 and cosine-shaped spin/charge 


density profile, as observed in numerical simulations H27H (de¬ 
tails to be given elsewhere). We are convinced that our funda¬ 
mental result shall find applications far beyond the treatment 
of boundary driven quantum master equation. For example, 
computer algebra (without a proof yet!) suggests an existence 
of an infinitely dimensional intertwiner (R-matrix) between a 
pair of auxiliary spaces, implying exact commutativity [13] of 
a two-parameter family of non-Hermitian transfer operators 

[Q(d, to), a(A', to’)] = 0, VT, A', to, to' e C. (28) 

We note that our novel Lax operator provides an appealing 
factorisation of our previous result IU6I1 - derived by a more 
brut-force approach - in the special case of zero spectral pa¬ 
rameter A — 0. We note that (A„i\(A, w)U =( O =o matches with 
the generator of Yangian symmetry of the Hubbard model 
[ 28|| truncated to a finite open chain of n sites. However, 
the relevance and facility of other conservation laws derived 
from C2(A, to), say for establishing rigorous bounds on trans¬ 
port coefficients |29ll, remain an exciting problem to study 
in future. We stress that our attempts to link our novel con¬ 
cepts to the 4-dim. Lax matrix constructed by Shastry 0], ot¬ 
to quantum- or super-symmetries of the related Hubbard-like 
models within the framework of integrable N — 4 super Yang- 
Mills theory [3(1 31 ], failed so far. It is thus in our opinion an 
urgent question to establish whether these results imply exis¬ 
tence of novel quantum symmetries of the ID Hubbard model. 
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